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Examples of nonsingular cosmological models are presented on the basis of exact solutions to multidimensional 
gravity equations. These examples involve pure imaginary scalar fields, or, in other terms, "phantom" fields with an 
unusual sign of the kinetic term in the Lagrangian. We show that, with such fields, hyperbolic nonsingular models 
with a cosmological bounce (unlike spherical and spatially flat models) emerge without special relations among the 
integration constants, i.e., without fine tuning. In such models, the extra-dimension scale factors as well as scalar 
fields evolve smoothly between different finite asymptotic values. Examples of theories which create phantom scalar 
, fields are obtained from string-inspired multidimensional field models and from theories of gravity in integrable Weyl 

■ space-times. 
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Introduction 



The recent years have been marked with a lot of violent events in both observational and theoretical cosmology. 
The discovery of an accelerated expansion of the Universe pQ has been one of the most important empirical findings. 
" . There followed a flood of theoretical works trying to interpret and to explain this acceleration, see, e.g., The 
majority of such constructions have a common feature: they involve various kinds of scalar fields. In many cases 
these are so-called phantom scalars, having a "wrong" sign of the kinetic term in their Lagrangians. 

It is of interest that such scalar fields are able to suggest a solution to one more long-standing problem of theoretical 
. £h ! cosmology, namely, the initial singularity problem. Inclusion of such fields makes it possible to circumvent the well- 
known singularity theorems and to prevent the formation of a cosmological singularity, keeping the curvature at 
H ' sub-Planckian scales. 



In these notes we will discuss in some detail this mechanism of avoiding a singularity in isotropic cosmological 
models at the level of classical field theory. Its efficiency is clear from the simplest example of a time-dependent 
scalar field in general-relativistic isotropic cosmology, described in Sec. 2. We begin with a discussion of the late-time 
behaviour of various models with accelerated expansion and show that a deceleration parameter qo < — I may be 
obtained with phantom scalar fields without creating a final singularity ("big rip") related to an infinite growth of 
the scale factor at finite time. It is then shown that an initial singularity can also be avoided with the aid of such 
fields: it turns out that a regular minimum of the scale factor is only possible in hyperbolic models with a phantom 
scalar field. 

A similar mechanism works as well in more complex cases to be discussed in the further sections. Sec. 3 gives 
examples of multidimensional models emerging in the field limit of some topical string-based unification theories 
(see ^ |S] and reviews on string cosmology [H]). Phantom type fields take place in such models with dimensions 
D > 11. We will not touch upon the whole wealth of exact solutions obtained in such models and restrict ourselves 
to simple models with a single dilatonic field, a single axionic type antisymmetric form and two scale factors — the 
"external" one, a(t), and the "internal" one, b(t). We shall see that nonsingular solutions in which, as t — > oo, 
the scale factor a(t) grows while b(t) and the dilaton <p(t) tend to finite limits, exist in closed and spatially flat 
cosmologies (Kq = 0, +1) with some special values of the integration constants only, i.e., require "fine tuning", 
whereas in hyperbolic models with Lobachevsky 3-geometry (Kq = — 1) they emerge generically, without any fine 
tuning. 

In Sec. 4, similar inferences are obtained for cosmological models of a certain class of gravitation theories with 
Weyl non-metricity which are, in many observational predictions, equivalent to scalar-tensor theories (STT) of gravity 
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and can also contain phantom scalar fields [S] . Their main difference from STT is a geometric interpretation of 
the scalar fields which makes their possible phantom character more natural than for usual, material scalar fields. 
Sec. 5 contains some concluding remarks. 



2. 4D scalar field cosmologies at late and early times 

To illustrate the general properties of scalar field driven cosmologies, let us consider the standard 4D Friedmann- 
Robertson- Walker (FRW) metric 

ds 2 = -dt 2 + a 2 (t)ds 2 , dag = 1 ^ + r 2 (d9 2 + sin 2 flrfy 2 ), (1) 

where a(t) is the scale factor and ds 2 ) is the metric of a 3-space of constant curvature Kq = 0, ±1 . The corresponding 
Einstein-Friedmann equations are 



a 



2 (a' 2 + K ) = xp, — (2aa" + a' 2 + K () ) = -xp, (2) 



where x — 8irG is the gravitational constant, the prime denotes d/dt, p = — Tq is the total density of matter and 
p = Tl = T% = T3 its pressure. The matter can be of any origin, but the symmetry of its stress-energy tensor 
T^j = diag(— p, p, p, p) is determined by the choice of the metric (JIJ. 

The most important kinematic observational parameters characterizing the Universe evolution are the Hubble 
parameter H(t) := a 1 /a and the deceleration parameter q(t) := —aa"/a' 2 . (The parameter q was introduced when 
it was believed that the expansion of the Universe was decelerating; its negative values correspond to an accelerating 
Universe.) The current observational estimates of these two parameters considerably vary from one paper to another 
but can be more or less reliably taken as 

ff «0.71±4 Jg-, q «-1±0.4, (3) 

s • Mpc 

where the subscript "0" refers to the present epoch. 

For p y^z 0, one can always write p — wp where w is in general time-dependent. The simplest models are obtained 
for w — const (the so-called barotropic matter). In this case the conservation law V a T^ = leads to the relation 

p = const -a- 3{w+1 \ (4) 

giving the well-known laws p ~ et~ 3 for dust (w = 0) and p ~ a~ 4 for disordered radiation (w = 1/3). Accelerated 
expansion (q < 0, however, requires a negative pressure. It is easily seen that, for w < —1/3 (which is needed for 
obtaining q < 0) and a large scale factor a, the term with Kq in J2J) is negligible as compared with p and p, so 
that the scale factor behaviour at late times does not depend on the spatial curvature. The t dependence of the scale 
factor is then described as follows: 

a) if —1/3 > w > —1 , so that the dominant energy condition holds, the expansion may be called power-law inflation: 

^ i a/[3(i»+i)] i q = _ 1 + l( w + !)>-!; (5) 

b) if w = —1 , which corresponds to a positive cosmological constant, p = const > 0, we obtain exponential inflation: 

a~e Ht , H = const, q = -1; (6) 



c) if w < — 1, the matter may be called exotic, and we obtain hyper- inflation ending with a singularity related to a 
blowing-up scale factor: 

o~(i,-t)- 2 /l 3 l ,0+1 l], 9 = -l-f|«; + l|<-l, (7) 

where i* is the singular time instant. 

In case c), matter behaves exotically indeed: its density grows as the volume grows, and all this ends with a "big 
rip" , a singularity at finite physical time, where both a and p grow infinitely. 
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Such a sad future of our Universe may be, however, avoided even if the present value of w is smaller than -1 but 
if w = p/p is time-dependent. Indeed, suppose that matter (or its dominating part) is represented by a scalar field 
(j) with the Lagrangian 

L s = -±£.9^,^ - V(<t>), (8) 
where e = ±1 and V(<fi) is a potential. For <j> = 4>(t) we have 

P=W + V, p=l £ ^-V, w=p/p = -i + —^— . ( 9 ) 

Thus a normal scalar field (e = +1) with a positive potential V gives w > — 1 wheras a phantom scalar field 
(e = —1) with V > leads to w < — 1. However, if at large t the scalar field tends rapidly enough to a minimum of 
the potential, U m i n > 0, then </>' — > and w — ► —1 as t — > oo; Vi n i n behaves as an effective cosmological constant, 
and accordingly we obtain the de Sitter asymptotic © ■ 

Let us now discuss a possible nonsingular behaviour of the early Universe, such that the scale factor a(t) undergoes 
a small but nonzero regular minimum a m j n at some instant t = to, so that a 1 (to) = 0, a" (to) > 0. Eqs. then give 

signp = K , p + 3p<0. (10) 

Suppose again that the model dynamics is dominated by a scalar field ||HJ| and that near a m i n the potential 
is negligible compared to the kinetic term, \V\ <C \4>' 2 \- Then, for <fi = 4>(t) , we obtain p p « £<t> 12 , and it 
is straightforward to find that the conditions 1)10(1 hold only in case Kq = —1, e = —1. Thus, among massless 
scalar fields, only the "phantom" one can lead to a bouncing isotropic cosmology, and only in hyperbolic models with 
Lobachevsky 3-geometry. 

One can note that the kinetic term, whose effective equation of state is p = p, evolves like a~ 6 [see (@J] whereas 
the potential term behaves qualitatively as a cosmological constant, therefore its neglection at small a(t) is justified. 



3. A multidimensional p-brane model 
3.1. General features 

Consider the action of D -dimensional gravity interacting with a dilatonic scalar field $ and antisymmetric forms F s , 
F p , which account for contributions from both the Neveu-Schwarz — Neveu-Schwarz (NS-NS) and Ramond-Ramond 
(RR) sectors: 



2 

D 



S, = / d u x^-g\ $ 



* 2 t-^i nJ s 



where 1Z is the scalar curvature, g = \ det gMN \ > (<9<I>) 2 = g MN dm&dN® , M, N = 0, . . . , D — 1 , u is a (Brans-Dicke 
type) coupling constant, n s and n r are the ranks of antisymmetric forms belonging, respectively, to the NS-NS and 
RR sectors of the effective action; for each n-form, F% — F n> M 1 ...M n ,Fjf 1 " Mn . 

The action i|ll|) is written in the so-called Jordan conformal frame where the field $ is nonminimally coupled 
to gravity. This form is actually obtained in the weak field limit of many underlying theories as the framework 
describing the motion of fundamental objects, therefore we will interpret the metric gMN appearing in (|ll|l as the 
physical metric. Thus, if the fundamental objects are strings, one has in any dimension lu — — 1 , while in cases where 
such objects are p-branes, one finds @] 

(D-l)(p-l)-(p+lf 

(D-2)(p-l)-(p + l)3' 1 ' 

where p is the brane dimension and D is the space-time dimension. The NS-NS sector of string theory predicts a Kalb- 
Ramond type field with n s = 3; the type II A superstring effective action contains RR terms with n r = 2, 4, while 
type IIB predicts n r = 3, 5. The action (|llfl may also represent the bosonic sectors of theories like 11-dimensional 
supergravity (where the dilaton is absent, and there is a 4-form gauge field), or 10-dimensional supergravity (there 
is a dilaton and a 3-form gauge field), or 12-dimensional "field theory of F-theory" [TJ]|, admitting the bosonic sector 
of 11-dimensional supergravity as a truncation. The model |10) contains a dilaton and two F -forms of ranks 4 and 
5; it admits electric 2- and 3-branes and magnetic 5- and 6-branes. 
The standard transformation 

9MN = <f- V{D - 2) 9MN (13) 
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leads to a theory reformulated in the Einstein conformal frame, more convenient for solving the field equations: 



(14) 



where all quantities are written in terms of the Einstein- frame metric gMN j ffE = | det Qmn I i f° r the scalar field we 
have denoted 

D — 1 ( D-V 



D 



while the coupling constants X s and A r are 
n. - 1 



A., 



wxC.D-2) 



(NS-NS sector); 



Ay 



77^ = sign lu 



2n r - L> 
2wi(D-2) 



D 



(RR sector). 



(15) 



(16) 



The sign factor i] u distinguishes "normal" theories {r) u = such that the kinetic term of the (p field in l|14(l 
has the normal sign corresponding to positive energy, from anomalous theories where this sign is "wrong" [rf u = — 1). 
The factor -q^ is thus quite similar to e used in Sec. 2. It should be noted that many theories with D > 11 involve 
f] w = —1. According to (JTSJ, 

(£>-2)(p-l)- 



^ = (£> - 2) 



'- ' (p+T)' ' J- (17) 

Evidently, under the condition (D — 2)(p — 1) > (p + l) 2 we have 77^ = — 1 . For p = 2, 5 this happens when £> > 1 1 , 
and for p = 3, 4 when £> > 10. 

The following table gives the values of u> and ij^/ujf for some particular space-time and brane dimensions. 



D 


P 








P 






any 


1 


-1 


L>-2 


12 


2 


-2 


-10/9 


10 


2 





8/9 


12 


3 


-3/2 


-5/2 


10 


3 


00 





12 


4 


-8/5 


-2 


10 


4 


2 


8/25 


12 


5 


-2 


-10/9 


10 


5 





8/9 


12 


6 


-6 


-10/49 


10 


6 


-4/9 


72/49 


12 


7 


1/2 


5/8 


11 


2 


CO 





12 


8 


-4/11 


110/81 


11 


3 


-2 


-9/8 


14 


2 


-4/3 


-4 


11 


4 


-5/2 


-18/25 


14 


6 


-16/11 


-132/49 


11 


5 


00 





26 


3 


-17/16 


-48 


11 


6 


1/4 


36/49 


26 


4 


-50/47 


-1128/25 



Some comments are in order. First, the well-known result to = —1 for strings (p = 1) in any dimension is 
recovered. Second, one obtains u) = 00 for 2- and 5-branes in 11 dimensions, which conforms to the absence of a 
dilaton in 11D supergravity that predicts such branes. Third, in 12 dimensions one has r\ u = — 1 for p < 7, and such 
a theory [10] does contain a pure imaginary dilaton: the F -forms of ranks 4 and 5 are coupled to a dilaton field tp 
with the coupling constants A 2 = —1/10 and A2 = — Ai , respectively, while the product Xip is real. As is concluded 
in Ref. ^U], for D > 11 "imaginary couplings are exactly what is needed in order to make a consistent truncation 
to the fields of type IIB supergravity possible" . In our (equivalent) formulation, ip and A are real and the unusual 
nature of the coupling is reflected in the sign factor rj^ . Supersymmetric models with D = 14 are also discussed 
[111 112] . while D = 26 is the well-known dimension for bosonic strings. 



3.2. Solutions 

There is a diversity of exact solutions (discussed, in particular, in Refs. |13l 1141 H3j . see also references therein) for 
the action (|14|l without L m , in space-times with the metric ~g AB of the form 

n 

dsl = - c 2Q ("W + J2 e 2/3 ' (u) ds 2 (18) 

i=0 

where u is a time coordinate and dsf are it -independent metrics of internal di -dimensional factor spaces assumed to 
be Ricci-flat for i = 1, . . . , n whereas ds^ describes the "external" (observed) space of constant curvature Kq = 0, ±1 , 
corresponding to the three types of isotropic spaces. 
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We will be only interested here in cosmological solutions for a very simple special case: a single antisymmetric 
form Fu ] from the NS-NS or RR sector, having a single (up to permutations) nontrivial component i*i...<2 where 
the indices refer to the external space Mo with the metric ds^ , a single internal space Mi with the metric ds\ , so 
that in Ijl8(l £ = 0,1, and (p — <p(u) . Then the field equations are easily integrated. 

Let u be a harmonic time coordinate for the metric l)18[l. so that the coordinate condition is a — do (3° + d\[3 1 . 

The F-form is magnetic- type; the Maxwell-like equations due to (|14f) are satisfied trivially while the Bianchi 
identity dF = implies 

Fi... do = Qy/go, Q = const, (19) 

where go is the metric determinant corresponding to ds^ and Q is a charge, to be called the axionic charge since 
the only nonzero component of F can be represented in terms of a pseudoscalar axion field in do + 1 dimensions. 
The remaining unknowns are (3° , (3 1 and ip . 

In the Einstein equations 1Z M — ^S^IZ = T^, written for the Einstein-frame metric (JTSJ, the stress-energy 
tensor Tj^ has the property + T* = (where z belongs to Mo ) , and the corresponding Einstein equation has 
the Liouville form a — (3° + K (d — l) 2 e 2a ~ 2/3 , whence 



do-i 



j3 - a = s(-K ,k,v) d = < 



e ku , Ko = 0, ke E; 

k^ 1 cosh ku, Ko = 1, k > 0; 

k~ 1 smhku, K = -1, fe> 0; (20) 

U, Ko = —1, k = 0; 

sinfett, ifo = — 1, k < 0, 



where fc is an integration constant and one more constant is suppressed by a proper choice of the origin of u . Eq. H20fl 
can be used to express j3° in terms of (3 = /3 1 . 

It is helpful to consider the remaining unknowns as a vector x A = ((3 1 , ip) in the 2-dimensional target space V 
with the metric 

couo -p*:). :)• <-) 

The equations of motion then take the form 



•• A 

X = -T} 2 



,Q 2 Y A e 2y (22) 

with the function y(w) = + A<p, representable as a scalar product of x A and the constant vector f in ¥: 

y(u) = Y A x A , Ya = (di, A), F A = (1/d, ^A). (24) 

Eq. IL'-it is a first integral of (1221 that follows from the ( u ) component of the Einstein equations. 
The simplest solution corresponds to Q = (scalar vacuum): 

(3 =cu + c, ip^c^u + c^, (25) 

where c 1 , c 1 , c v and are integration constants. Due to (|23J) . the constants c" 4 = (c 1 , c v ) are related by 

c A c A = dd 1 (c 1 ) 2 + i lu cl = —^—K. (26) 



If Q 7^ , Eqs. 1221) combine to yield an easily solvable (Liouville) equation for y(u) : 

y + VF Q 2 Y 2 e 2 y = 0, Y 2 = Y A Y A = d x /d + VoJ X 2 . (27) 

This is a special integrable case of the equations considered, e.g., in Refs. ^3E1E3- Assuming 2 Y 2 > 0, Eq. lt271) 

gives 

e -3/(Vj =h- 1 \Q\Ycosh[h(u + u 1 )} (28) 



2 Even for r/^ = — 1 one has Y 2 > for fields from the NS-NS sector in any dimension and for fields from the RR sector if D < 17. 
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where Y = \Y 2 \ 1 / 2 , h > and u\ are integration constants. The unknowns x A are expressed in terms of y as 
follows: 

Y A 

(29) 



x" = —r-yfu) + c u + c A 
Y z 



c A Y A 



where the constants c A — (c 1 , c v and c A — (c 1 , c ) satisfy the orthogonality relations 

0, c A Y A = 0. 

Finally, the constraint leads to one more relation among the constants: 

do 



Y 2 



c A c 



do-1 



-K. 



(30) 



(31) 



3.3. Analysis of cosmological models 

In what follows, we put do = 3, so that d\ = D — 4, and identify, term by term, the Jordan-frame metric ds 2 
obtained in the above notations iJTSJl, l(T%j) . 



ds 2 



exp 



' 2i r 



2S(-K ,k,u) 



-du 2 



AS 2 {-K 0l k,u) 



ds 2 



^ds\ 



where the function S(. 1 ., .) is defined in (|20|l . with the familiar form of the metric 

ds 2 = -dt 2 + a 2 (t)dsl + b 2 (t)dsl, 



(32) 



(33) 



so that a(t) and b(t) are the external and internal scale factors and t is the cosmic time. 

To select nonsingular models, let us use the Kretschmann scalar JC — RmnpqR RINP ® , which is in our case a sum 
(with positive coefficients) of squares of all Riemann tensor components Rmn P ® ■ Thus as long as K, is finite, all 
algebraic curvature invariants of this metric are finite as well. For the metric (|33|l with d = 3 one has (the primes 
denote d/dt): 



K = A 



3 - +d i r- 



., , ,'a'V 



2-i 



Kn 



(34) 



By (|34f) , K, — > oo and hence the space-time is singular when a — > , a — * oo , b or b — > oo at finite proper 
time t. Accordingly, our interest will be in the asymptotic behaviour of the solutions at both ends of the range 
K« = («min, « m a) of the time coordinate u, defined as the range where both a 2 and b 2 in (|33|) are regular and 
positive. (Note that, as t — * ±oo, a singularity does not occur when b(t) — > 0, or a — * in case Ko = 0.) At any 
u G M u all the relevant functions are manifestly finite and analytical. The boundary values u max and it m i n may be 
finite or infinite; a finite value of u max or u m i n coincides with a zero of the function (|20|l . 

Among regular solutions, of utmost interest are those in which a(t) grows while b(t) tends to a finite constant 
value as t — > oo . Any asymptotic may on equal grounds refer to the evolution beginning or end due to the time- 
reversal invariance of the field equations. We will for certainty speak of expansion or inflation, bearing in mind that 
the same asymptotic may mean contraction (deflation). 

Let us now enumerate the possible kinds of asymptotics. 



Type I: 



±oo , 



dr 



M-2k)\u\ du 2 



(35) 



with k > and the constants 4 an B, depending on the parameters of the solution. An asymptotic of interest for 
a(t) takes place if A > 2k > 0: 



(i) A > 2k: 

(ii) A = 2k: 



t — > oo , a ~ 

t A/{A-2k) 

(power-law inflation) ; 
t ~ |u| — > oo, a ~ e kt (exponential inflation). 



A reformulation for k < is evident. The scale factor b(t) tends to a finite limit if B = 0, i.e., under a special 
condition on the model parameters (fine tuning). 

Type la: a modification of type I when k = 0, so that at u — > oo 



df 



.-3 „A 



e Au du\ 



-1 „A 



e Au du 2 , 



(36) 
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If A > , we have, as desired, t — ► oo and a — ► oo ; the expansion may be called "slow inflation" since it is only 
slightly quicker than linear: the derivative da/dt ~ it, which behaves somewhat like Inf. If A < 0, then a — ► at 
finite t (singularity). As for b{t) , one may repeat what was said in case I. 

Type II: u — > 0, where the function (|20|) tends to zero, so that S(—Ko,k,u) ~ u, while other quantities involved 
are finite. In this case 

dt 2 ~l/u 3 , a 2 ~l/u, b 2 ^ const > 0. (37) 

According to JHZJl, f — > ±oo, a(t) ~ |t| (linear expansion or contraction), whereas both b(t) and ip(t) tend to finite 
limits since they do not depend on S(—Ko, k, u) . 

The dilaton tp in all cases behaves like In b(t) , but, in general, with another constant B in each particular solution. 

This exhausts the possible kinds of asymptotics for Y 2 > 0. Solutions with Y 2 < 0, which can emerge when 
r] u = — 1, may have other asymptotics, but they are of lesser interest. 



Scalar-vacuum cosmologies 

The scalar- vacuum models <|32[) . (|25[) depend on two input constants, D (or di = D — 4, or d = (D — 2)/2) and 
a; (or wi) and three integration constants k, c , c v related by (|26|) : two more constants, c 1 and c^, only shift the 
scales in Mi and along the ip axis and do not affect the qualitative behaviour of the models. 

Closed models, Kq = +1 . In this case in 1-iL'l) S — k/ coshfcu, k > 0, hence the solution has two type I asymptotics 
at u — > ±oo , with k > and the following constants A — A± : 



A± = -k T 



d^ + 



(38) 



so that at least at one of the asymptotics A < whence a — > at finite i , a singularity. The behaviour of 6(f) is 
also singular. 

Spatially fiat models, Kq — 0. One has simply 

a 2 (t)=e Au , dt^ e iA - 2k ^ 2 du, (39) 

where A = —c v / {du)\) — die 1 — k, k e R, and again 6 2 (f) = e B ", S = const. Thus each of the scale factors is either 
constant, or evolves between zero and infinity, and a = occurs at finite t. 

Hyperbolic models, Kq = — 1. If k > [note that, when rj^ = 1, there is necessarily fc > due to l|2fi[l]. one has 
in (|32(l S — k^ 1 sinhfcu. Hence the model evolves between a type I asymptotic at u — > oo, with A coinciding with 
A + in Eq. Ij38(l , and type II at u = . Since type II is regular, a necessary condition for having a nonsingular model 
is A>2k. 

To find out if and when it happens for r)uj = +1, it is convenient to introduce, instead of the two constants c 1 
and c v connected by H26(l . an "angle" 6> such that 



~ cl = \ I^TT k cos8, -Ctp = J^-k sm6. (40) 

V 2dd\ V 2 

The condition A> 2k will be realized for a certain choice of the integration constants if A + given by l|38|l has, as a 
function of 9 , a maximum no smaller than 2k . An inspection shows that it happens if 

W? < l/[d(6d - di)] = 1/[(U - 1)(U - 2)]. (41) 

This is the only example of a nonsingular (bouncing) vacuum model with -q^ = +1. 

In case k > 0, r/^ = — 1 , a choice of c v and c 1 subject to (|26|l such that A > 2k is easily made for any u>t . 
For r] u = — 1, fc = 0, the model evolves between type la and II asymptotics, where at the la end (u — * oo) 

A = -die 1 - c v /(duji), B = 2c x — c v /(duj\). (42) 

The necessary condition for regularity, A > , is satisfied for proper c\ and c v which can be chosen without problems. 

In case = — 1, k < 0, the function (|20|l is simply sin \k\u, and the model has two type II asymptotics 
at adjacent zeros of S , say, u = and it = 7r/|fc| . This model is automatically nonsingular for any further choice of 
integration constants. 
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We conclude that among vacuum models only some hyperbolic ones can be nonsingular. For r/^ = +1 in such a 
case a(t) evolves from linear decrease to inflation, or from deflation to linear growth. Only in the latter case both 
b(t) and ip tend to finite limits as t — * oo without any fine tuning. 

For r\ u = — 1 there is a model interpolating between two asymptotics of the latter kind. Thus, as t changes 
from — oo to +oo, a(t) bounces from linear decrease to linear increase (generically with a different slope) whereas 
b(t) and (p(t) smoothly change from one finite value to another. The latter model exists for generic values of the 
integration constants. 



Cosmologies with an axionic charge 

The solution contains, in addition to the input parameters D, oj and A, three independent essential integration 
constants: the "scale parameter" k, the charge Q and also h and c v connected by ((31 (1 : the constant c 1 is excluded 
by the first relation (|30|l 

die 1 + Xc v = (43) 

so that the quantity c a ca, appearing in (|31|l . is expressed as c a ca = Vuj(d/di)c^Y 2 . The fourth constant, the "shift 
parameter" u\, as well as c 1 and c^, connected by JSOJl, are qualitatively inessential. 
Let us begin with "normal" models, r\ u = +1. The solution 1)29(1 has the form 



1 (u) = -^y(u) + c 1 u + c 1 , y{u) = -^y{u) + c v u + c^. (44) 



where the function y(u) is fiven by J35J). The condition (|3f \ leads to k > and strongly restricts the possible model 
behaviour. Thus, it can be shown [3] that for all K one of the asymptotics belongs to type I with the constants 



dYA 1+ ui J ' dun. 



A = -k + — [ di + ] + c v ( A 



The necessary condition for regularity A > 2k may be fulfilled for small u>\ , satisfying the condition l|41|) . just as in 
the vacuum case. The constant B may then have any sign, and only a special choice of the ratio c v /h (fine tuning) 
can lead to B — , providing a finite limit of b(t) . 

The second asymptotic depends on Ko . For closed and flat models it is again type I, and a nonsingular behaviour 
is again acheived by fine tuning. For Kq = — 1, the second asymptotic belongs to type II and is always regular. Thus 
particular models with an axionic charge may be even regular for r) u = 1 ) and, in addition, may be inflationary as 
t — > oo. 

The "anomalous" models with r) u = —I, just as in the vacuum case, are more diverse due to arbitrariness in the 
sign of k. For models Kq = 0, +1 , as well as for hyperbolic ones with k > 0, the restriction 1)41(1 is no more valid, 
but the type I and la asymptotics are again only nonsingular for special values of the parameters. 

Lastly, the models Kq — — 1, rj^ = — 1, /c<0,as their vacuum counterparts, interpolate between two type II 
asymptotics and have the same qualitative features. 



4. Vacuum multidimensional models with integrable Weyl geometry 

Let us discuss multidimensional cosmological models assuming that space-time possesses D -dimensional Weyl geom- 
etry characterized by the metric gMN and the connection 

r A . c = f A . c -±(<j B 5£,+a c 6 A -g B c<T A ) (46) 

where Tg C are the Christoffel symbols for the metric gAB , u is a scalar field and a a — 9a<J . The gravitational field 
is described by the tensor gAB and the scalar a , as in scalar-tensor theories (STT). Just as in STT, the gravitational 
Lagrangian may contain different invariant combinations of gAB and a . Restricting ourselves to Lagrangians which 
are linear in the curvature and quadratic in a a , we can write: 

L = f(a)K - h{a)a A a A ~ 2A(<r) + L m (47) 

where 1Z is the Weyl scalar curvature, obtained from the connection 1(46(1 . /, h and A are arbitrary functions, and 
L m is the non-gravitational matter Lagrangian. 



Nonsingular multidimensional cosmologies with Lobachevsky spatial sections 







The field equation are simplified if one expresses the Weyl curvature 1Z in terms of the Riemannian curvature R , 
corresponding to the metric g A B '■ 

TZ = R-(D -lpff -^(D-l){D-2)cr A a A (48) 

( R and □ are formed from the Riemannian connection Tg C ) and passes to the Einstein conformal picture with the 
aid of the transformation gMN — f~ 2 ^ D ~ 2 ^9~MN- Omitting a total divergence, we arrive at the following form of 
the Lagrangian: 

L = R-F(a)g AB a A a B + f- D /( D - 2 \~2A(a) + L m ] (49) 
where R is the Riemannian scalar curvature for the metric ~g A B > A a = dA/da and 

(50) 



D-2" 



Consider vacuum (L m — 0) cosmological models with the metric (|18fl . assuming A = 0. One can easily find 
that the substitution a i— ► ip, such that dtp /da = y/\F(a)\ , leads the action with the Lagrangian 149(1 to a form 
coinciding with l(14jl without an F-form and with the sign factor ry w replaced with signF(cr). Therefore, in the 
Einstein picture, the vacuum cosmologies of the gravitation theory with Weyl integrable space-time are entirely 
identical to the scalar-vacuum models from Sec. 3 (both types of models are equivalent to those of multidimensional 
general relativity with a massless, minimally coupled scalar field). 

A difference can appear after a transition to Jordan's picture since the conformal factors $" 2 /( £>_2 ) in Sec. 3 and 
f- 2 /( D - 2 ) i n the present section are, in general, different. However, if one supposes that the function f(a) is finite 
and smooth in the whole essential range of a , then the qualitative properties of the models (as regards their regular 
or singular behaviour) coincide in the Einstein and Jordan pictures. Accordingly, in both pictures, the conclusion 
that hyperbolic models in the presence of a phantom scalar field contain a class of nonsingular bouncing models, 
preserves its generality. Such models do not require any fine tuning, and in all of them both the scalar field and the 
extra-dimension scale factors change in finite limits. 

Stability of the qualitative features of these models with respect to addition of other kinds of matter has been 
confirmed by a numerical study of Weyl cosmologies with one more scalar field representing ordinary matter . 



5. Concluding remarks 

Phantom scalar fields are rather widely discussed as one of the dark energy candidates, able to explain the present 
accelerated expansion of the Universe, see Sec. 2 and, in more detail, e.g., |16| . Such fields, if any, may also dominate 
in the early Universe, at small values of the scale factor a(t), above all, due to maximum stiffness (p = p) of the 
equation of state of a massless scalar field, so that their energy density p grows with falling a(t) more rapidly than 
for other kinds of matter: e.g., in 4D FRW models one has p ~ a~ 6 . For the same reason, massive scalar fields, or 
fields with potentials, should actually behave as massless ones at small a(t) . Therefore the above conclusion that an 
initial cosmological singularity may be avoided without fine tuning in hyperbolic models (Kq = — 1) due to phantom 
scalar fields seems to be rather general. 

This conclusion proves to be even more important in multidimensional cosmologiies. If extra dimensions, being 
an inevitable ingredient in modern unification theories, are considered dynamically, the singularity problem becomes 
even more involved since, in addition to the usual cosmological scale factor, the extra dimensions can collapse or 
blow up, leading to a curvature singularity. However, in open models of the type discussed here, the external scale 
factor a(t) dynamically differs from other variables. Formally, this circumstance was described here by the last 
line in Eq. (|20[1 : a 2 (t) ~ 1/S(—Kq, fc,it) tends to infinity at two finite values of the harmonic coordinate u, which 
correspond to infinite physical time, whereas all other variables, such as scalar fields and internal scale factors, remain 
finite at these values of u. 

The above simple models certainly do not pretend to describe the whole evolution, but only try to guess the 
qualitative features of the bouncing process near the maximum density state. The well-known features of standard 
cosmology: inflation, nucleosynthesis, particle creation etc. may follow at later times, but the mechanism described 
here seems to automatically provide stable compactification of the extra-dimension scale factors (if any) and constant 
values of scalars which may be related to coupling constants in unification theories, such as the dilaton in string 
theory. Cosmological and astrophysical problems related to stable compactification are discussed in Ref. ^7], see also 
references therein. 

I would also like to mention that such an exotic matter as a phantom scalar field, violating all standard energy 
conditions, in case it is concentrated in comparatively small regions of space, is precisely what is needed to create 
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wormholes — and hence maybe also time machines |18| . There are many exact wormhole solutions involving phantom 
scalars (those in Refs. ^HlGHl are probably the earliest). If such fields do exist in nature, whatever be their origin, 
they are, in principle, a ready construction material for wormholes. 

Observations are known to yield the total cosmological density factor Qq smaller or close to unity; meanwhile, 
the presently popular spatially flat cosmologies, most convenient for various calculations, require the precise equalty 
Q,o = 1 , actually a sort of fine tuning. It is much more probable that the real Universe at least slightly violates this 
special requirement, leading to Kq = — 1 if Qq < 1. 
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